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The inclusion rating method by statis- 
ti cs of extreme values (IRMSE) using 
\Zarea of inclusions as the size parame- 
ter enables one to discriminate between 
current super-clean steels. Moreover, 
IR MSE enables one to predict the size 
(y/area m,,) of maximum inclusions con- 
tained in domains larger lhan the in- 
spection domain. The statistical 
distribution of ^area„„ can be used 
for the quality control of materials and 



for the prediction of a scatter band of 
fatigue strength. Practical procedures of 
inclusion rating and prediction of a 
scatter band of fatigue strength arc 
shown. 
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1. Introduction 

With the increase in cleanliness of steels, con- 
ventional inclusion rating methods are no longer as 
useful as before, because conventional inclusion 
rating methods cannot determine the cleanliness of 
new clean steels. Although the cleanliness of steels 
has been markedly improved in the last two 
decades, the fatigue strength of recent clean high 
strength steels cannot attain the ideal value ex- 
pected from their high static strength. Nonmetallic 
inclusions are predominantly the cause of lower fa- 
tigue strength even for such clean high strength 
steels. Thus, in order to predict the fatigue strength 
behavior and to evaluate quality, we need a new 
inclusion rating method relevant to recent super- 
clean steels. The inclusion rating method based on 
statistics of extreme [1] is most relevant for this 
purpose. In the following, we call this method In- 
clusion Rating Method by Statistics of Extreme 
(IRMSE). 

In this study, we shall first show that if we choose 
an appropriate size parameter for inclusions, the 



size of inclusions obey the statistics of extreme 
value theory. The appropriate size parameter is the 
square root of projected area of the maximum in- 
clusion co ntain ed in a standard inspection area or 
volume, y/area mxx . Second, we predict the size of 
the maximum inclusion which may be contained in 
a larger area or volume than the standard inspec- 
ti on ar ea and, lastly, we use the size parameter, 
^/arcana, to predict the scatter band of fatigue 
strength of hard steels. 

The merits of IRMSE, in comparison with con- 
ventional methods, are (1) to distinctly discrimi- 
nate the cleanliness of recent super-clean steels, 
and (2) to predict the size of larger inclusions con- 
tained in a domain larger than the inspection do- 
main. This method is useful for quality control of 
materials and for improvement of the steel making 
processes. It also enables one to predict the scatter 
of the fatigue strength of a large number of mass 
production products. 
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2. Nonmetallic Inclusions as a Fatigue 
Fracture Origin 

Figure 1 shows an example of the nonmetallic 
inclusion which was observed at fatigue origin of a 
bearing steel under a rotating bending fatigue test. 
If this inclusion did not exist in this specimen, the 
fatigue strength of this specimen should have been 
higher than the applied stress, <r a = 1078 MPa. 
Since the size and location of nonmetallic inclu- 
sions scatter randomly, the fatigue strength of high 
strength steels naturally scatters. Although there 
has been a firm opinion that the chemical composi- 
tion and shape of nonmetallic inclusions substan- 
tially influences the fatigue limit, Murakami et al. 
[2 5] have shown the incorrectness of the conven- 
tional opinion by their detailed experiments and 
analyses, and reported distinct experimental evi- 
dence that the size of inclusions (defined by 
y/area) is the most crucial geometrical parameter. 
Tt is empirically known that the intrinsic fatigue 
strength of steels is determined by the hardness 
(Hv) of its microstructure. For steels with // v <400, 
nonmetallic inclusions contained in current com- 
mercial steels are not detrimental and we have the 
following empirical formula 



o*w — 1>6 Hv 



(1) 




where or K is the fatigue limit (MPa) and Hv is the 
Vickers hardness (kgf/mm 2 ). However, for steels 
with // v >400, the effect of inclusions reveals itself 
and the intrinsic or ideal fatigue limit given by Eq. 
(1) cannot be attai ned. The fatigue strength de- 
pends on the size (y/area) and location of the fatal 
inclusion and H v of the matrix. Murakami et al.'s 
[6-9] fatigue limit prediction equations are classi- 
fied into three categories depending on the loca- 
tion of fatal inclusions (see Fig. 2): 

Fatigue limit for a surface inclusion [Fig. 2(a)] 



<r« - 1.43 (H v + 120)/(vWa ) 



1/5 



(2) 



Fatigue limit for an inclusion in touch with free 
surface [Fig. 2(b)] 

tr w = 1.41 (H v + 120)/(Vowfl) !/ft . (3) 

Fatigue limit for an internal inclusion [Fig. 2(c)] 

ow = 1.56 (Hv + 120)/( ySreu )'*, (4) 



where the units are <r M : MPa, \Jarea : [im, and H v : 
kgf/mm 2 . 



Free surface 
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Fig. 2. Various locations of inclusions causing fatigue fracture. 

Since for a constant value of area, an inclusion is 
most detrimental when it exists just in touch with 
the free surface of a specimen, we can use Eq. ( 3) 
in combination with the maximum size y/area max 
obtained by IRMSE to predict the lower bound 
(ow) of scattered fatigue strength of many speci- 
mens or machine elements. 



Vickers hardness Hv = 745 kgf/mm : 

Applied stress at surface o-„ =t078 MPa 

Cycles to failure M=7.94 x lO* 

Square root of projection area of inclusion \/areu = 9.8 u,m 

Distance from surface h = 158 u,m 

Applied stress at inclusion o-' = 1034 MPa 

Chemical composition of inclusion: AI-Mn-S-0 

Fig. 1. A typical example of inclusion observed at the center of 
fatigue fracture origin [super-clean bearing steel, SUJ2(N)]. 



3, Inclusion Rating of Various High 
Strength Steels by Statistics of 
Extreme 

Figure 3 explains the practical procedure to im- 
plement the inclusion rating by statistics of extreme 
values. The details of this method are reported in 
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/ areamKj(^u m) 




M.999S 



So(mm 2 ) 



Fig. 3. Practical procedure of the inclusion rating by statistics of 
extreme values. 

Murakami et al.'s papers [3,5,9-1 1]. The procedure 
is briefly explained in the following, see [11]. 

(1) A section perpendicular to the maximum 
principal stress is cut from the specimen. After pol- 
ishing with a n°2000 emery paper, the test surface 
is mirror-finished with buff. 

(2) A standard inspection area So (mm 2 ) is fixed. 
Generally, it is advisable to take a microscope pic- 
ture for an area approximately equivalent to So. In 
the area So, the inclusion of maximum size is se- 
lected. Then, the square root of the projected area 
vWoma of this selected inclusion is calculated. 
This operation is repeated n times (in n areas 5 U ) 
(see Fig. 2). 

(3) The values of y/area^aj are classified, start- 
ing from the smallest, and indexed: (with; = 1..ji ). 
We then have the following relation: 



y/area mia ^ y/area^i* 



■\/area B 



The cumulative distribution function F t and the 
reduced variates y, are then calculated from the 
equations. 

F, =j X 100/(n + 1) 

y ; = -In[-!n07(/H-l)]. 

(4) The data are then plotted on probability pa- 
pe r. T he point j has an abscissa coordinate of 
y/areanvxj while the ordinate axis represents either 
Fj oTy h An example of the curve is shown in Fig. 4. 

Figure 4 shows the inclusion ratings by IRMSE 
for two kinds of super-clean bearing steels, 
SUJ2(N) and SUJ2 (H). The total oxygen con- 
tained in these steels is 8 ppm for SUJ2(N) and 
5 ppm for SUJ2(H). This kind of information 
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Fig. 4. Cumulative frequency of the extreme values of inclusions 
[Super-clean bearing steels, SUJ2(N) and SUJ2(H)]. 

enables one to discriminate quantitatively the dif- 
ference among the cleanliness levels of the same 
kind of materials produced by different companies 
or produced by a company at different periods. 
Thus, this information will be useful for the quality 
control of materials and the improvement of the 
steel making process. 

It is not a pri ori e vident to what extent the ex- 
treme values \/area m3I . of inclusions contained in 
various steels follow extreme statistics value. How- 
ever, Murakami et al. [3,5-11] have shown many 
examples of measurements which obey the statis- 
tics of extreme value theory. Uemura and Mu- 
rakami [12] carried out a three-dimensional 
numerical simulation to find t he sta tistical distribu- 
tion of the extreme values \/area m3I of inclusions 
which were distributed in a constant volume with 
the size (D) distribution of the type, 



<KD)=-exp(--), 



where m is the mean value, and they confirmed the 
validity of IRMSE (Fig. 5). In addition, they indi- 
cated the quantitative difference between two-di- 
mensional and three-dimensional measurements, 
though the difference virtually vanishes with in- 
creasing inspection domains. 
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(a-1) Inclusion distribution on 
inspection section of the 
material with m=l pxn. 
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(b-1) Inclusion distribution on 
inspection section of the 
material with m=2 Mm. 
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(C-l) Inclusion distribution on 
inspection section of the 
material with m=3 um. 
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(b-2) Statistical distribution of 
true(3-D) and apparent(2-D) 
extreme values of inclusions 
for material with m=2 |xm 
(So=0.25 mm2» number of 
inspections = 40). 
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(c-2) Statistical distribution of 

true(3-D) and apparent(2-D) 
extreme values of inclusions 
for material with m=3 um 
(So=0.25 mm 2 , number of 
inspections = 40). 



Fig. 5. Numerical simulation of the inclusion rating by statistics of extreme values on the materi- 

1 D 
als with the inclusion size distribution of the type <£(D) = — exp ( ). 
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4. Application to Prediction of Scatter 
Band of Fatigue Strength 

Figure 6 illustrates the shape and dimension of a 
tension-compression fatigue specimen [13]. The 
material used is tool steel, SKH51. The chemical 
composition is shown in Table 1. Table 2 shows the 
mechanical properties. 
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Fig. 7. Statistical distribution of the extreme values, the maxi- 
mum size of inclusion at the center of fracture origin (Tool 
steel, SKH51). 



Fig. 6. Shape and dimension of tension-compression fatigue 
specimen (mm) (Tool steel, SKH51). 

Figu re 7 shows the extreme value distribution of 
\Jarea of the inclusions found at the fracture origin 
of 34 specimens. The data in Fig. 7 are the extreme 
values obtained by the fatigue test but not by the 
two-dimensional metallographic method described 
in Sec. 3. Figure 8 indicates the location of these 
inclusions on the fracture surface. If the tension- 
compression fatigue test is not performed correctly, 
that is, specimens are subject to a bending moment 
due to a bad alignment or the curving of the speci- 
men axis, nonmetallic inclusions existing near the 

Table 1. Chemical composition in wt% of material (Tool Steel, 
SKH51) 



IjO 



c 


Si 


Mn 


P 


S 


Cr 


W 


Mo 


0.81 


0.31 


0.29 


0.018 


0.002 


3.92 


6.10 


4.85 


V 


Co 


Cu 


Ca 


AI 


Mg 


o 




1.81 


0.46 


0.07 


0.004 


0.035 


0.0005 


0.0018 





100 



II 



o <?m- MPa, Specdia.«9mm 
a "m- 784MPa,Spec.dia.-9mai 
a Ojn"-784MPa, Spcc.dia.-6mm 



50 C 






0.1 0.2 0.3 0,« 0.5 0.6 0.7 0.8 0.9 1.0 

Distance from surface/Radius of specimen 

Fig. 8. Relationship between the size {y/orta) and location of 
inclusions at the center of fracture origin (Tool steel, SKH51). 

free surface are likely to appear as the fracture 
origin on the fracture surface [14]. In such a case, 
unusually low fatigue strength is likely to be ob- 
tained. Since the fatigue fracture origins shown in 
Fig. 8 are distributed randomly on the section of 
specimen, these data may be valid for the statistical 
analysis. However, it should be noted that when 



Table 2. Mechanical properties of quenched and tempered test material (Toot steel, 
SKH51) 
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the surface inclusions became the fracture origins, 
the data were not plotted on Fig. 7, because such 
inclusions are more detrimental than an inclusion 
having the same size and existing internally and ac- 
cordingly they may be a little smaller than the exact 
maximum inclusion. 

In the case of the data of Fig. 7, the volume of 
the test part of one specimen (Fig. 6) corresponds 
to one inspection domain and there are 34 extreme 
values in Fig. 7. Therefore, Fig. 7 can be used for 
predicting the expected maximum size of the inclu- 
sion which may be contained in more specimens 
than those used in fatig ue tests. For example, an 
inclusion having y/area m ^=\3^ u,m is expected to 

be contained in 100 specimens (iv" = 100). 

Combining this \/aream^ ( = 138 |xm) and Eq. (3), 
the lower bound (a- wf ) of fatigue strength of 100 
specimens can be predicted. 

Figure 9 compares the scatter observed in exper- 
iments and the predicted lower bound o-„c of the 
scatter band. The prediction is in good agreement 
with experiments. The prediction of the lower 
bound of fatigue strength explained above can be 
used for the quality control of machine elements 
which are produced by mass-production and can- 
not be tested individually. 

The data as shown in Fig. 7 offer us reliable in- 
formation on inclusions expected to be contained 
in other specimens. However, obtaining the data 
shown in Fig. 7 requires preparation of many pre- 
cise specimens and time consuming fatigue tests. 
To avoid this inconvenience, the author has pro- 
posed an alternative two-dimensional method as 
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Fig. 9. Comparison between the experimental results and the 
lower bound of fatigue strength which was predicted on the ba- 
sis of Eq. (3) and the maximum size of inclusion (Tool steel, 
SKH51). 



explained in Sec. 3. A sufficient number (TV) of in- 
spection domains (insp ection areas) necessary to 
predict reliably \/area mM for more specimens or 
larger areas should depend on the materials to be 
inspected and on the inspection area S a observed 
by the image processor combined with an optical 
microscope. From the author's experience, it is rec- 
ommended that N be larger than 40 for So = 0.031 
mm 2 . 

Several Japanese industries have already put the 
method proposed in this study in practice [15]. 



5. Conclusions 

(1) If we define the size of nonmetallic inclu- 
sions contained in commercial stee ls by the square 
root of th e pr ojected area, y/area, the maximum 
values, y/area m „, in a definite inspection domain 
obey the statistics of extreme value theory. 

(2) The inclusion rating method by t he sta tistics 
of extreme values (IRMSE) based on y/area^^ can 
be used for a new inclusion rating method. IRMSE 
enables one to discriminate distinctly between re- 
cent super-clean steels, while conventional inclu- 
sion rating methods are no longer valid as the 
method to evaluate the cleanliness of new clean 
steels. 

(3) IRMSE is useful not only for a relative eval- 
uation of materials but also for the prediction of 
the expected maximum size of inclusions to be con- 
tained in a domain la rger than the inspection do- 
main. The value of y/area^ can be used with the 
fatigue strength prediction equation to predict a 
scatter band of fatigue strength of high strength 
steels. 
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